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We model the 3K-phase of Sr2Ru04 with Ru-metal inclusion as interface state with locally 
enhanced transition temperatures. The resulting 3K-phase must have a different pairing sym- 
metry than the bulk phase of Sr2Ru04, because the symmetry at the interface is lower than in 
the bulk. It is invariant under time reversal and a second transition, in general, above the onset 
of bulk superconductivity is expected where time reversal symmetry is broken. The nucleation 
of the 3K-phase exhibits a "capillary effect" which can lead to frustration phenomena for the 
superconducting states on different Ru-inclusions. Furthermore, the phase structure of the pair 
wave function gives rise to zero-energy quasiparticle states which would be visible in quasiparti- 
cle tunneling spectra. Additional characteristic properties are associated with the upper critical 
field Hc2- The 3K-phase has a weaker anisotropy of Hc2 between the inplane and z-axis orien- 
tation than the bulk superconducting phase. This is connected with the more isotropic nature 
Ru-metal which yields a stronger orbital depairing effect for the inplane magnetic field than in 
the strongly layered Sr2Ru04. An anomalous temperature dependence for the z-axis critical 
field is found due to the coupling of the magnetic field to the order parameter texture at the 
interface. Various other experiments are discussed and new measurements are suggested. 

KEYWORDS: Sr2Ru04, p-wave superconductor, broken time-reversal symmetry, frustration effects, Andreev 
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§1. Introduction 

In recent years the quasi-two-dimensional metal 
Sr2Ru04 has advanced to one of the most intensely stud- 
ied transition metal oxides displaying unconventional su- 
perconductivity.^' There is strong evidence for spin- 
triplet pairing with broken time reversal symmetry, a 
pairing state with the basic form d(k) = z{kxiiky) (chi- 
ral p-wave state). '^"^^ This superconducting phase shows 
a variety of unusual properties among which the recently 
discovered "3-Kelvin" (3K) phase is one of the most puz- 
zling findings. While the transition to the bulk super- 
conducting state occurs around 1.5K, in samples with a 
large excess-Ru concentration a precursors to the super- 
conducting transition appears at temperatures as high 
as 3K.^'^^ This phase shows the features of an inhomo- 
geneous superconducting phase. The detailed material 
analysis revealed that the excess-Ru does not distribute 
uniformly, but forms small inclusions of micrometer-size. 
Thus, the 3K-phase is likely connected with the phase 
separation of Ru-metal and Sr2Ru04 in this sample. The 
bulk superconductivity of metallic Ru has a supercon- 
ducting transition temperature of 0.5 K only with a con- 
ventional (s-wave) pairing state. This leads us to the 
assumption that the boundaries between the two mate- 
rials provides the environment for the local nucleation of 
superconductivity at a higher temperature (Fig. 1). 

It is not our aim to discuss here in detail the origin for 
the locally enhanced transition temperature, since the 
microscopic theory of superconductivity in Sr2Ru04 is 
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still unclear. Nevertheless, we would like to comment on 
one important aspect which could be connected with the 
enhanced superconductivity. Various experiments have 
shown that a particular soft optical modes for local lat- 
tice distortion is associated with the inplane rotation of 
the RuOg-octahedra, the Ss-mode at the Brillouin zone 
boundary of the phonon spectrum, which leads to a slight 
volume reduction.*^ This rotation affects one of the three 
electron bands in particular, the 7-band which originates 
from the 4d-i2g-orbital with (ia;j,-symmetry.'^'^) While the 
dispersion for this orbital occurs via 7r-hybridization be- 
tween the Ru-d-orbital and the 0-p-orbitals in the undis- 
torted case, the rotation introduces an additional a- 
hybridization of opposite sign for the da;j,-orbital. A lat- 
tice distortion of this kind would diminish the dispersion 
of the 7-band and enhance the electron density of states, 
because of a Van Hove singularity near the Fermi level. 
It is likely that the Ru-inclusions in Sr2Ru04 lead to in- 
ternal stress that is released by local static distortions in 
the vicinity of the interfaces, most likely connected with 
the Ea-mode {Uujs^ ~ 140K). A crude estimate from 
neutron scattering data yields a length scale of order 50 
- lOOA over which this rotational distortion extend away 
from the interface. The increased density of states would 
lead to a locally enhanced independent of the micro- 
scopic mechanism. The additionally enhanced ferromag- 
netic spin fluctuations may support the spin-triplet pair- 
ing instability too.^"^ 

In this paper we would like to investigate a number of 
properties of the inhomogeneous 3K-phase from a phe- 
nomenological point of view. A generalized Ginzburg- 
Landau (GL) formulation is most suitable for this pur- 
pose, since we will discuss an inhomogeneous supercon- 
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ducting phase. The basic assumption is that the inter- 
face region has an enhanced transition temperature. We 
briefly review the basic conclusions of our theory. The 
first and most important fact is that the superconducting 
state nucleated at the interface has a different symmetry 
than the bulk superconducting phase. The interface su- 
perconducting state is time reversal symmetry conserv- 
ing. Consequently, there is a further second order phase 
transition where this symmetry is broken. Wc will show 
that this transition occurs in general above the onset of 
bulk superconductivity. Naturally the nuclcation of su- 
perconductivity at the interface is inhomogcneous and 
does not lead to a uniform phase transition. This is also 
true for the second transition. Even if all interfaces be- 
tween Ru-inclusions and Sr2Ru04 are locally equivalent, 
their geometry and mutual arrangement would lead to 
a spread on nucleation temperatures due to "capillary 
effects". 

The 3K superconducting state at the interface corre- 
sponds to an odd parity state, a p-wave state with a 
pair wave function that has a node parallel to the nor- 
mal vector with a positive and negative lobes parallel to 
the interface (Fig. 1). This phase structure is responsi- 
ble for a peculiar change of the quasiparticle spectrum 
due to Andreev reflection, i.e. the accumulation of An- 
dreev bound states at zero-energy which may be observed 
by quasiparticle tunneling. The phase structure together 
with the topology of the interfaces can generate a frustra- 
tion of the order parameter phase which can be released 
by introducing spontaneous orbital currents. Note that 
also the second transition to a time reversal symmetry 
breaking phase would reduce the frustration and intro- 
duce spontaneous currents. 

Finally the upper critical field Hc2 shows special prop- 
erties in case of fields parallel to the interface. In this case 
the orbital depairing is reduced due to the essentially two 
dimensional nature of the condensate. In addition, the 
anisotropy of Hc2 between fields parallel to basal plane 
of the strongly layered Sr2Ru04 and the 2;-axis direc- 
tion is reduced compared with the critical field for the 
bulk phase. This is due to the more isotropic nature 
of the Ru-inclusions. A considerable fraction of the su- 
perconducting condensate of the 3K-phase resides in the 
Ru-part. An additional very striking feature occurs for 
magnetic field parallel to the z-axis of Sr2Ru04. While 
in zero field only the p-wave component with its node 
perpendicular to interface appears, the magnetic field 
drives also the other component whose node is parallel 
to the interface. The energy gain occurs via the coupling 
of the field to the current induced by the order parame- 
ter texture of both components. This coupling leads to 
a peculiar enhancement and temperature dependence of 
Hc2 in agreement with experiment.^' 

§2. Phenomenological description 

The following analysis is based on the generalized GL 
theory as the most efficient way to describe basic prop- 
erties of an inhomogcneous superconducting state. The 
bulk pairing state of Sr2Ru04 has the symmetry of a 
chiral p-wave state, represented by d(k) = Aoz{kx ± 
iky). This requires a two-component order parameter 



V = iVx^ily) with d(k) = z(t7 • k), which belongs to 

the two-dimensional representation E-^ of the tetragonal 
point group D^h and describes the leading instability in 
Sr2Ru04. Note that recent flux distribution measure- 
ments in the mixed state have suggested the presence of 
two order parameter components. On the other hand, 
Ru is a conventional s-wave superconductor with a tran- 
sition temperature around 0.5 K. The superconductivity 
in Sr2Ru04 penetrates the Ru-metal due to the proxim- 
ity effect in two ways. First, proximity leads naturally to 
a spin-triplet pairing amplitude in Ru, although its crit- 
ical temperature there may be extremely small. Second, 
the spin-triplet superconducting state in Sr2Ru04 can 
induce the s-wave component. The corresponding cou- 
pling is, however, probably weaker than for the triplet 
channel because at the interface the triplet and singlet 
spin wave function have to be connected, by means of 
spin-orbit scattering. We will ignore the s-wave compo- 
nent in Ru and will briefly comment later only. 

2. 1 Ginzburg- Landau free energy 

The GL free energy for the two-component order pa- 
rameter T] has the well-known form, 

+K2{\Dy7]^f + \D^Vyf) + {K3{D^ri^)*{Dyriy) 

+Ki{{Dyr]^)*{D^r]y)+c.c.} 

+K5{\D,r],\'' + |£>.%|') + (V X A)V87r] 

(1) 

where the coefficients are different in the two subsystems 
which we label by indices S and R for Sr2Ru04 and Ru, 
respectively. The gradient terms contain the gauge in- 
variant spatial derivatives D = V -|- i{2e/hc)A with A 
denoting the vector potential. The coefficients Ki^j_ deter- 
mine the coherence length of the superconducting order 
parameter (/z = S and R). Since Sr2Ru04 has a lay- 
ered structure, the coherence length along the z-axis is 
short, K^s ^ Kis,K2S, ■■■■ On the other hand, Ru is 
more isotropic so that K^n has a magnitude similar to 
the other coefficients. In a weak-co\ipling approach as- 
suming cylindrical or spherical Fermi surface shapes we 
obtain the relation Kifj_/3 = = K^fi = ^5^4^. The 
second order coefficient af^{T) changes sign at the bare 
bulk transition temperature T^^. For our discussion of 
the qualitative properties of the interface superconduc- 
tivity it is sufficient to assume linear temperature de- 
pendence an(T) = a^iT — Ten) "^'^^^ ^cS ~ 1.5K and 
TcR = 0. 

For simplicity, we consider a single homogeneous pla- 
nar interface with a normal vector parallel to the x-axis 
(n = (100)) where — oo < x < belongs to the Ru- 
metal, and < .X < oo to Sr2Ru04. In this geometry 
the problem reduces to a one dimensional with the spa- 
tial direction along the a;-axis. In the following we will 
always ignore the z-direction assiiming homogeneity in 
this direction for simplicity. On the Sr2Ru04-side a thin 
layer of thickness d with enhanced transition tempera- 
ture is introduced. The boundary conditions, in general, 
involve reflection and transmission of Cooper pairs at 
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the interface. We make the simplifying assumption of 

complete transparency, which corresponds to a continu- 
ous order parameter. We will, however, comment on the 
more general case below. We may describe the thin layer 
at the interface by a (5-function in the free energy, if d 
is much smaller than the coherence lengths along the x- 
axis (Fig. 1). Thus, the interface part of the free energy 
has the form, 



(2) 



where a = daa{T — Tea) with Tea > TcS- This approxi- 
mation is sufheient to discuss most important qualitative 
features of the 3K-phase. 

SroRuQ, 




Fig. 1. Interface between Sr2Ru04 and a Ru-inclusion. The in- 
terface has a layer of thickness d of ciiiianced transition tempera- 
ture where a p-wave state nucleates whose wavefunction has the 
lobes parallel and the nodes perpendicular to the interface. 



2.2 Instability conditions of a planar interface 

We will first investigate the conditions of the nucle- 
ation of superconductivity at the planar interface. We 
consider a temperature range Tcs < T < Tea where this 
local onset of superconductivity is supposed to occur. 
The instability condition is obtained by searching for the 
solution of the linearized GL equations which are given 
by 



Kip.dlvx - a^Vx = 0, 



(3) 
(4) 



for x 7^ (/z = S and R for x > and x < 0, respec- 
tively). At X = the solutions have to be continuous 
and satisfy the following boundary conditions, 

Klsda;T]x\x=0+ - Kijid^-nxlx^o- - (^'nx\x=o = 0, (5) 
K2sdx'ny\x=0+ - ^2r9x%U=0_ - cr%|x=0 = 0. (6) 

For temperatures above Tcs the order parameter is 

largest at the interface and decays exponentially on both 
sides. For the ?7a;-component the solution is 

Vx{x) = ??a;oe""=/^^^ for x > (7) 
=??xoe"/^^=^ forx<0 (8) 



with = Kifjt/a^. The analogous solution exists for 
Tjy with = K2^j./a,^. Then Eq.(5) and (6) lead to the 
instability equations. 



\/ K2sa.s ^ 



2RaR + O" 



(7 = 

0, 



(9) 
(10) 



for rjx and ijy, respectively. Because Ki^ > > 
(Ci/i > Cl/i) the nucleation will occur for the y- 
componcnt. To obtain an instability temperature T* 
above T^s requires that 



C2R{Tcs)aR{Tcs) < daa{Tca - T 



cS) 



(11) 



This simple relation gives a good insight on the basic 
problem of the nucleation of local superconductivity at 
the interface. Obviously, the larger d and the higher 
Tea, the higher is T*. However, the presence of the 
normal-metal Ru tends to suppresses superconductivity, 
in particular if the coherence length ^2R increases. Note 
that if both sides of the interface were Sr2Ru04, we al- 
ways would find local superconductivity at a tempera- 
ture higher than Tcs- Furthermore, the suppression by 
the Ru-metal is weaker, if the interface were less trans- 
parent leading to a discontinuity of the order parame- 
ter. Reduced transparency of the interface would not 
change the dominance of the r^y-component at this in- 
terface. The reason is that the order parameter compo- 
nent n • is suppressed by a reflecting interface. In any 
case the degeneracy of rjx and rjy of the bulk region is 
lifted at the interface, since it corresponds to a region of 
effectively reduced symmetry. The dominance of the 
r]y component is also plausible from a microscopic point 
of view. The opening of a gap for momentum direc- 
tions along the interface leads to the gain of condensa- 
tion energy for quasiparticles with momenta parallel to 
the interface. These are the quasiparticles spending the 
longest time in the interface region. 

We assume from now on that the condition Eq.(ll) 
is satisfied and the resulting transition temperature is 
T* « 3K corresponding to the S/C-phase. Note that the 
transition temperature does not depend on the orienta- 
tion of the normal vector as long as it lies in the basal 
plane. The nucleating order parameter t] is perpendicu- 
lar to the normal vector, i.e., n x r/. The situation does 
not change much if a small ^;-axis components of the nor- 
mal vector is introduced and the properties of the local 
superconducting state is determined by its inplane com- 
ponents. 

Because the 3K-phase does not break time reversal 
symmetry in contrast to the bulk phase below 1.5K, a 
further second order phase transition has to occur at a 
temperature T2 which is strictly larger than the bulk 
Tg. This is a different concept from conventional in- 
homogeneous superconductors, where superconductivity 
nucleates locally on small "islands" which with lowering 
temperature increase their overlap and, finally, form a 
bulk superconducting phase without further symmetry 
breaking (apart from the percolation transition). The 
temperature Tj* depends on the coupling between the 
order parameter components described by the fourth or- 
der terms in the free energy. The two pairing components 
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may "attract" or "repel" each other, determined by the 

sign of the parameter b = 2bi — &2 + ^3 > repul- 
sion and 6 < attraction) . The weak coupling approach 
leads to 6 > 0, which suppresses the appearance of the 
second component. Hence Tj* may be rather close to 
and difficult to distinguish experimentally from the bulk 
transition. 

The discussion of the second transition and the form 
of the order parameter in the 3K-phase is complex as it 
involves also vector potential due to the presence of spon- 
tancoTis ciirrents in the time reversal symmetry-breaking 
phase. Thus we turn to the numerical solution of the 
complete set of GL equations including the complete set 
of GL equations for the order parameter and the vec- 
tor potential. We use the parameters as — ckr = eta- 
and Tea = 2.8TcS. Except for K^fj^ we choose all coeffi- 
cients Kjf^ to be the same in Ru as well as in Sr2Ru04: 
Kifj^/a^^Tcs = iK^^^/a^Tcs = 1 (i = 2,3,4), which cor- 
responds to the basic zero-temperature coherence length 
^0 ~ 1- The coefficient for the z-axis gradient is small 
in Sr2Ru04, K^^/a^Tcs = 10"^, while we take, for sim- 
plicity, JTsR = K2K, since Ru is more isotropic. The 
fourth order coefficient are also chosen independent of 
X, 2bif_i ~ 362^ = —363^ = OAaij,. We choose the inter- 
face layer thickness d = 1, larger than it probably is in 
reality, for illustrative purpose and numerical stability. 
With these parameters we obtain a nucleation temper- 
ature of T* « 2r,;s and T2 » 1.33rcs which is rather 
high due to the large value of d . At this transition the 
rix component becomes finite in addition to rjy and has 
the relative phase of ±7r/2. In Fig. 2 we show the shape 
of the two order parameter components in the three dif- 
ferent temperature regimes. Here a) represents the gen- 
uine 3K-phase with vanishing ry2:-component (time rever- 
sal symmetry conserving phase), b) is the intermediate 
time reversal symmetry breaking state with both com- 
ponents finite. Finally, c) is the bulk superconducting 
phase of Sr2Ru04. Note, that both cases b) and c) pos- 
sess a complex order parameter texture at the interface 
which will be important for the magnetic properties dis- 
cussed in the next section. 



2.3 Capillary effect 

The interface instability for the superconducting state 
bears some resemblance with the wetting of a surface. 
Similar to the wetting phenomena of liquids we find here 
capillary effects. Modulations of interfaces on length 
scales comparable to the coherence length can enhance 
the nucleation temperature. While this capillary effect 
is rather simple for conventional superconductors, there 
are complications in the case of an unconventional super- 
conductor. In particular, there can be frustration effects 
due to the internal structure of the pair wave function, 
as we will see in the next section. Here we would like to 
consider first the rather simple situation of an interface 
that is not flat but has a weak modulation. 

The enhancement of the transition temperature on a 
spatial modulated interface can be most simply inter- 
preted by the effective increase of surface area. This 
corresponds also to an enhanced mutual overlap of the 




X 

Fig. 2. Spatial dependence for the two order parameter compo- 
nents dashed line and \riy\ solid line): a) 3 K -phase at 
T = 1.75Tcs > Tg; b) intermediate time reversal symmetry 
breaking phase at T = 1.25TcS < with both components 
which have a relative phase ±7r/2 two-fold degeneracy); c) 
bulk superconducting phase d = z{kx ± iky) and a texture at 
the interface at T = 0.95TcS. 



order parameters nucleating at different points of the in- 
terface or on different inclusions. Let us assume that 
the interface is only slightly modulated, described by 
xo{y) = lsm{2-jTy/L') where d <^ I <^ L' ^ ^i^. Then, 
by the most simple variational approximation we replace 
the interface term by 

:F, = J d'r\rj\Ul+l^^co.l^^^^%{x-xo{y)). 

(12) 

The important modification appears via the new inter- 
face metric, which accounts for the fact that the inter- 
face is wider or "denser". This leads to an enhanced 
instability temperature, even if we approximate the spa- 
tial dependence of the order parameter by a form like in 
Eq.(7,8). Using this variation approach the instability 
equation for r]y changes to 

\/K2sas + VK2Ran + aR = 0, (13) 
where the factor R is 

(14) 

Obviously, R is always larger than 1 and leads to an 
effective enhancement of o". We have ignored in our vari- 
ational approach the spatial dependence of the order pa- 
rameter along the interfacc\ Including this aspect lead to 
even further enhancement of the nucleation temperature. 

Another aspect of the capillary effect is the mutual in- 
fluence of interfaces which we would like to consider on 
the example of two parallel interfaces. To be specific we 
assume a thin Ru-metal slab of thickness L sandwiched 
between Sr2Ru04, again with normal vector parallel to 
the a;-axis. Both interfaces have the same properties. 
Then it is easy to derive the instability equations taking 
the boundary conditions into account. There are two 
combinations of the nucleating order parameters on the 
two interfaces: a "bonding" and "antibonding" configu- 
ration which is even or odd, respectively, under reflection 
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at the center of the slab. For the instability equations of 
the r/y-component wc obtain, 



^ K2sas + tanh 



tinction between magnetic properties of the 3K-phase at 
T* and the bulk phase. 



26r 



VK2RaK + a = (15) 



for the bonding and 



V K2sas + coth ( ) V K2RaR + a = 



(16) 



for the antibonding combination. Obviously, we recover 
in both cases the original instability equation Eq.(lO), 
if we separate the two interfaces far apart, L ^2R- 
If, however, L ^ ^2R then the second (positive) term 
is diminished (enhanced) and a higher (lower) transition 
temperature T* results in the case of bonding (antibond- 
ing). Naturally, the same kind of capillary effect occurs 
also in the inverse situation where a Sr2Ru04 slab is 
surrounded by Ru. The bonding combination of the two 
interface states just corresponds to the adjustment of the 
phases of two superconducting islands. The antibonding 
combination is equivalent to a phase difference of tt and is 
the energetically least favorable case. We will sec in the 
next section that this aspect is important for frustration 
effects of the 3K phase. 

Wc can conclude that the onset of the 3K phase is 
rather inhomogeneous on the interfaces as well, because 
the capillary effects have a strong influence on the nucle- 
ation of the superconducting order parameter. 

§3. Physical consequences 

3.1 Spontaneous interface magnetism 

While the onset of superconductivity is visible in the 
reduction of electrical resistance, the significant features 
of second transition at T2 are less easy to detect. Bulk 
superconducting double transitions are often observed 
through specific heat anomalies. Here the only signifi- 
cant anomaly, however, is seen at the onset of bulk su- 
perconductivity at Tc = 1.5K. Magnetic properties may, 
on the other hand, allow us to observe the onset of the 
superconducting phase at T2 that violates time reversal 
symmetry. The solution of the full GL equations show 
that the below T| spontaneous supercurrents occur par- 
allel to the interface. In Fig. 3 wc show the magnetic field 
and current distribution at the interface for the case cor- 
responding to Fig. 2b. There are two currents along the 
interface flowing in opposite direction as a result of the 
texture of the two order parameter components at the 
interface. Screening currents are weak due to absence of 
superconductivity in the bulk. Thus, the field resulting 
from the currents is strongly peaked on the interface and 
can even generate an overall finite magnetic fiux. Since 
the flux is finite, it could be observed, in principle, by 
high-resolution magnetic microscopes, such as scanning 
Hall probes or a SQUID microscopes. Another sensitive 
probe for local magnetic field distributions is zero-field 
muon spin relaxation, which is expected to show an in- 
crease of the depolarization rate below T2. However, if 
T2 is very close to Tc (as is very likely the case) then the 
signal will be obscured by the bulk superconducting state 
which has previously been observed by /xSR.^^ The muon 
spin relaxation experiment would in any case give a dis- 
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Fig. 3. Spontaneous magnetic field (solid line) and supercur- 
rent jy (dashed line) distribution for the state given in Fig.2b) 
at T = 1.25TcS. 



3.2 Frustration effects in the inhomogeneous state 

We have seen in the previous section on the capillary 
effect that the lowest energy configuration of order pa- 
rameters nucleated on different interfaces naturally cor- 
responds to equal order parameter phases. Let us now 
see how the superconducting state would arrange on sev- 
eral inclusions close enough to each other that their order 
parameters substantially overlap. If the order parame- 
ter were conventional, it would be easy to adjust all the 
order parameter phases on the different inclusion to be 
equal. This is not the case, if the order parameter has 
p-wave symmetry as is illustrated in Fig. 4. The three 
Ru-inclusions depicted in Fig. 4 carry a superconducting 
state which on every point of the interface correspond 
to a p-wavc state with the momentiim direction aligned 
with the interface, i.e. the gap node lies always parallel 
to the normal vector. The interface regions on different 
inclusions close to each other behave like junctions or 
weak links connecting a network of the superconducting 
islands. 

The direction of the positive lobe of the pair wave func- 
tion is indicated by the arrows in Fig. 4 where we assume 
that the order parameter phase is constant on each in- 
clusion. Parallel (antiparallel) arrows on neighboring in- 
terfaces correspond to bonding (antibonding) configura- 
tions of the overlapping order parameters, or equal phase 
(TT-shifted phases) . It is obvious that in the arrangement 
of the order parameter on each inclusion introduces a 
TT-phase shift between its weak links. Although we may 
change the phases on each island, it is impossible to ad- 
just them so that all links have a zero phase difference for 
loops consisting of an odd number of inclusions (Fig. 4 
shows this for the case of three inclusions) . The situation 
is in many respects similar to the case described for gran- 
ular (i-wave superconductors.^^' We may identify 
loops which contain tt phase shifts and are consequently 
frustrated, since they cannot adjust the order parameter 
phase to minimize the energy of all weak links simulta- 
neously. The frustration can be released by introducing 
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phase gradients (supercurrents) which are energetically 

favorable, if the weak links grow strong enough. This 
would yield orbital magnetic moments. This effect re- 
sults from the competition between weak link energy and 
magnetic field energy (determined by geometry) as in 
usual frustrated 7r-loops. However, the frustration also is 
reduced with the time reversal symmetry breaking tran- 
sition at r = Both ways have similar impact by di- 
minishing frustration and lead to spontaneous currents. 
In the network of superconducting islands the transition 
to the time reversal symmetry breaking state is inhomo- 
geneous. 




Fig. 4. Configuration of the p-wave pair wave function on Ru- 
inclusions (sliadowcd regions) in a crosscction parallel to the 
a;-3/-plane. The arrow indicates the direction of positive wave 
function (see box). The dashed circles show the regions form- 
ing weak links between inclusions. The fact that the arrows of 
different inclusions in the weak link region are antiparallel cor- 
responds to a phase difference of tt. In the given configuration 
there is no way to adjust the phases of all three inclusions so 
that all weak links have minimal energy, i.e. a vanishing phase 
difference. Therefore this system is frustrated. 



The frustrated 3K-phase with spontaneous currents 
can exhibit enhanced absorption in the ac-susceptibility. 
Dissipative processes, such as phase slips, are associated 
with the hysteretic reversal of the spontaneous currents. 
Thus the absorption would be sensitive to the applica- 
tion of a small static external fields which would bias the 
spontaneous currents (a similar situation was observed 
in granular high-temperature superconductors^^' ^*''^''^) 
The presence of spontaneous orbital magnetic moments 
in the inhomogeneous time reversal symmetry violating 
phase above TcS can cause a characteristic non-linear 
magnetic response in small external fields, similar to 
the paramagnetic Meissner effect or Wohlleben effect in 
high-temperature superconductors.^^' ^^'^"^ Under field- 
cooling conditions the magnetic moments would be gen- 
erated in a polarized way yielding a paramagnetic con- 
tribution. Probably, the signal would be considerably 
weaker than in high-temperature superconductors. The 
inhomogeneous 3K-phase is, however, more suitable to 
measure the spontaneous magnetism than the bulk su- 
perconducting phase, where the macroscopically visible 
magnetism is only resulting from surface currents, while 
it is screened or compensated in the interior of the sys- 
tems. 



3.3 Quasipartide spectrum 

The presence of superconductivity modifies the quasi- 
partide spectrum in small enclosed normal-metal regions 
through the formation of so-called Andreev bound states. 
Andreev bound states are standing waves of an electron- 
hole superposition, e.g. in a normal metal region en- 
closed by a superconductor. In a quasiclassical picture 
the standing wave corresponds to an electron and hole 
travel on the same classical (ballistic) trajectory, but in 
opposite direction and are subject to Andreev reflection 
upon impact in the superconductor at both ends of the 
trajectory. The energy of such a standing wave depends 
on the phase of the gap functions at the boundaries of 
the classical trajectories and can be estimated by Bohr- 
Sommerfeld quantization. Energies much smaller than 
the gap of the superconductor are approximately given 
by, 

E{n,x)^^[{2n+l)n + x] (17) 

where vp is the Fermi velocity, L the length of the trajec- 
tory, n an integer and, most important, x is the difference 
between the phases of the gap functions at both ends of 
the ballistic trajectory (for simplicity we ignore here the 
effect of impurity scattering which can also be the origin 
low-energy states). For unconventional superconductors 
the phase difference x depends also on the direction of 
the electron (hole) momentum, because the gap func- 
tion, in general, has an anisotropic phase structure, i.e. 
a different phase for different directions of momenta on 
the Fermi surface. Thus, even in case of a homogeneous 
order parameter, the phase x appearing in the Andreev 
scattering process would no trivial, but revealing the in- 
ternal phase structure of the pairing state. For real order 
parameters we find x = which is most important, be- 
cause it leads to a zero-energy state independent oi vp 
and L. This phase difference appears, if the gap function 
has positive and negative sign depending on the momen- 
tum direction on the Fermi surface. From the quasiclas- 
sical point of view underlying Eq.(17) the geometry of 
the classical trajectory decides about the phase x = or 
TT. Since a large fraction of trajectories can have x = 
and all of them yield a zero-energy state, the zero-energy 
states (x = tt) lead to an enhanced density of quasiparti- 
de states at the Fermi surface. The presence of this kind 
of zero-energy bound state has been intensively investi- 
gated for surfaces of rf-wave stipercondcutors.^^'^^^ 

In our case a p-wavc pairing state nucleated at the 
interface generates a gap function which has regions of 
different sign on the Fermi surface. Zero-energy Andreev 
bound states are likely to occur within the Ru-inclusions 
as well as in Sr2Ru04 between Ru-inclusions that are 
sufficiently close to each other. We encounter here a 
complex form of an (imperfect) Andreev billiard. The 
Andreev reflections are not perfect. Due the finite width 
of the superconducting regions, the zero-energy level ac- 
quires some width. Nevertheless, there is a strongly en- 
hanced density of states at zero energy. Once time rever- 
sal symmetry is broken the phase differences x for tra- 
jectories deviate from and n. The originally enhanced 
density of states at zero-energy spreads over a larger en- 
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ergy region. (Note that the reduction of the density of 

states at the Fermi energy represents also a driving force 
for the time reversal symmetry breaking transition, as 
was suggested also for the surface states of a d-wave su- 
perconductor.^®' ^^)) 

One experimental indication for the enhanced density 
of states zero-energy states is the observation of so-called 
zero-bias anomaly, i.e. an increased quasiparticle tun- 
neling conductance at zero voltage. Recent measure- 
ments of quasiparticle tunneling conductance in c-axis 
facing break junctions of Sr2Ru04 with Ru-inclusions 
report the observation of zero-bias anomalies.^®-' In the 
3K-phasc a zero-bias anomaly in current-voltage charac- 
teristics develops gradually with decreasing temperature. 
This conductance peak deforms into a pronounced bell 
shape combined with a residual zero-bias anomaly, be- 
low the bulk superconducting transition. Assuming that 
these tunneling features reflect the quasiparticle density 
of states connected with the Andreev billiard, we can in- 
terpret them within the scenario of the nucleated p-wave 
state which turns into a time reversal symmetry breaking 
state at lower temperature. We would like to emphasize 
the fact that break junctions in samples without Ru- 
inclusions did not show any similar features above and 
below the onset of bulk superconductivity at 1.5 K. Note 
also that an analogous phenomenon, the deformation 
of the zero-bias anomaly, has been observed for tunnel- 
ing conductance into the [110]-surface of YBa2Cu307.^®-* 
Also in that case broken time reversal symmetry is most 
likely cause. ^^'^^^ 

It is clearly desirable to have more extensive exper- 
imental investigation of the quasiparticle spectrum in 
and around the Ru-inclusions via tunneling and related 
probes. ^"^"^^^ In particular, using a contact to a sin- 
gle inclusion and measuring the quasiparticle current- 
voltage characteristics through the Ru-metal inclusion 
into Sr2Ru04 could provide further valuable informa- 
tion on the 3K as well as bulk superconducting phase. 
A detailed analysis of the phenomena discussed in this 
section will be given elsewhere. 

§4. Upper critical field 

The 3K phase has an upper critical field which exceeds 
the bulk critical field considerably and shows a weaker 
anisotropy of in-plane versus out-of-plane critical field. 
The fact that superconductivity is confined in a small 
layer at the interface naturally leads to an reduction of 
orbital depairing, if the external magnetic field is applied 
parallel to the interface. This type of effect is known for 
superconductivity at other planar defects such as twin 
boundaries. ^^'^^^ In a magnetic field superconductivity 
nucleates first at interfaces which are parallel to the field. 
We will analyze this effect for the [100]-interface with 
magnetic fields B _L x-axis. Unfortunately, a complete 
analytical treatment of this problem is not possible even 
for the GL formulation. Therefore we will restrict our 
analytic discussion to the region of very small fields at 
the onset of the 3K-phase and use a variational approach. 
This allows us to illustrate a few basic features of the 
upper critical field. Then we will consider the behavior of 
Hc2 by numerical means for the model introduced above. 



For temperatures close to T* the critical field is small 
and we may consider it as a perturbation.^'*^ Thus, we 
will use the exponential form of the order parameter 
appearing in zero-field at the onset of the 3K-phase: 
Vxf^ = ?7rroexp(-|.T|/^i^) and 77^^^ = r]yoe-Kp{-\x\/(^2fi) 
with = Kifj_/a^. The magnetic field lies in the plane 
of the interface B = iJ(O,sin6',cos0) with the vector 
potential A — H{x — xo){0, — cos 6, sin 6). While in zero- 
field only the ?7y-component is nucleated at T*, in general 
both components can appear in a finite magnetic field. 
We use now the given form of the order parameter and 
calculate the free energy in a finite field by integrating 
over the spatial coordinates. In this first-order pcrtur- 
bative form the free energy per unit area of the interface 
is then up to second order in the the order parameter, 

l2 ^ , ^ , 7 



E 

M=S,R 



(7 

2 ' 4 



"2 ^?2^/w(^) 

+'i{v*xoVvo - Vxovlo)nHK2^£,2^ cos 9] 

(18) 

with 7 = 2e/nc and r = {Vs - l)/(\/3 + 1)^ where 
we used the relation = ■\/3^2^, keeping the weak- 
coupling relations among different Kif^ for cylindrical 
symmetry. To simplify the free energy (18) we have fixed 
xq to zero, since it plays a minor role for our discussion. 
The anisotropy parameters are 

fxf.iO)=K2^cos-'e 



Kr,^ sin^ 9, 
K^^, sin^ 9. 



(19) 



First we consider the case of the in-plane field {9 = 
7r/2). Here only the ryj^-component appears, which corre- 
sponds to the polar state with its nodes perpendicular to 
the field, which also appears in the bulk superconducting 
phase for in-plane fields. The zero of the corresponding 
coefficient in defines the instability, which leads to 



c2 



-d 



(20) 



where we used that Kr^s ^ Kifj_, Kr^u (i 7^ 5). In this case 
including a finite value of xq we find that xq = — ^2fl/2 
maximizes Hc2 (we omit here an explicit demonstration). 
Note that we can write Ep(C2A<a^(T) +cr(r)/2) = a*(T) 
which is proportional to (T — T*) close to T*. Con- 
sequently, we observe a square-root dependence Hc2 oc 
|T* — T|2.23,24) Pqj. (juj-vcd interfaces the regions tan- 
gential to the field tend to allow the nucleation at higher 
temperature. Simultaneously, we have, however, to in- 
clude capillary effects. Thus, the observation of a pure 
square-root temperature dependence would be masked 
in reality. 

For the field parallel to the z-axis {9 = 0) the rjx- 
component is involved too, despite the lower critical tem- 
perature in zero field. The optimal value of xq is very 
close to zero so that we fix xq = 0. With this simplifi- 
cation the instability equations involving the coupling of 
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Fig. 5. Upper critical fields for the 3K and bulk phase: The nu- 
merical solution of the GL equations give the upper critical fields 
of the 3K-phase: inplane (circles) and z-axis (diamonds). The 
corresponding bulk critical fields are given by doted lines. Note, 
the weaker anisotropy of the critical field of the 3K-phasc. The 
numerical data compare well with the experimental data on a 
qualitative as well as quantitative level.®' '^^ 



both components have the form, 



ijHA'rj^o + {<(T) + Ay}r]yo = 



(21) 



where 



Moreover, 



^"ImCIm and A' 



+ 



'^Ep^ZM^M- Moreover, <(T) = J:^,{^l^.a^{T) 
a(T)/2) is the effective second order coefficient for the 
77a:-component and changes sign at a temperature Tex 
which Hes between and T*. The critical field is ob- 
tained by searching the non-trivial solution of this equa- 
tion system. 



Ay a*(T) -4A'VA, 



(22) 



where we used \a*y{T)\ < a*(r) for T ^ T*. Thus 
we find again a basic square root behavior of Hc2(T) as 
in the case of inplane fields. However, in addition we 
see that the coupling of the two components yields an 
enhancement factor to i?c2 which becomes stronger as 
the temperature is lowered. (Note, however, that the 
present form is only valid for < a%iT) 4A'^/Aj, and 
in any case our variational approach is applicable only 
in a restricted region close to T*.) This enhancement is 
important, since it may modifies the overall form of the 
temperature dependence of Hc2- Before considering this 
point numerically let us, however, address the issue of 
Hc2 anisotropy. 

For the bulk superconducting phase of Sr2Ru04 the 
critical fields are given by 



Hi, 



-$oas(r) 



and H, 



-$oas(r) 



2S 



c2 



(23) 



27rVA'2sA'5S 

with T < Tcs and ignoring anisotropy of Hc2 for differ- 
ent inplane orientations, which is small. This leads to 
an anisotropy factor of order Hf.2/H^ « 2V^2s/v'-K'5S 
whose experimental value is about 12 for T close to 



Tcs in experiments {K2S > K5S ~ 0.03K2s).^^ On the 
other hand, the anisotropy for the 3K-phase may be less 
anisotropic with a factor 4 ~ 5,^'^^ since 



^2R^5R 



(24) 



gives a smaller ratio due to the fact that Ru-metal is 
basically isotropic with K^n ^ K2f_i and also the coher- 
ence lengths at T = T* are of the same order in Ru and 
Sr2Ru04. The reduction of anisotropy originates from 
the isotropy of the Ru-mctal which leads to a stronger 
coupling of the order parameter to the inplane field than 
for Sr2Ru04. 

The numerical evaluation of the upper critical field for 
the model used in Sect. 2. 2 illustrates the temperature 
dependence on qualitative level. The result are shown 
in Fig. 4 for the assumption that /fss = 0.03X23 and 
/sTsR = K2S ■ The circles and diamonds are the numerical 
results for the critical field inplane and along the 2;-axis, 
respectively. We observe indeed a weaker anisotropy for 
the critical field of the 3K-phase than for the bulk phase 
also indicate in Fig. 4. The initial temperature depen- 
dence has square root dependence as obtained analyti- 
cally. While the inplane critical field has a downward 
curvature in the plotted temperature range, the ^;-axis 
field is nearly linear. The analytical expression in Eq.(22) 
indicates that even the upward curvature would be pos- 
sible for the 2:-axis critical field. Our result is close to the 
experimental data qualitatively and to some extent quan- 
titatively, with the choice of parameters giving T* m 2Tc 
and the proper anisotropy of the bulk critical fields.^' 
The peculiar temperature dependence of the z-axis crit- 
ical field of the 3K-phase is also similar to that found in 
experiment.^' 

§5. The s-wave order parameter 

Ru metal becomes a conventional superconductor at 
a transition temperature of about 0.5 K. Thus, we may 
include an s-wave order parameter into our theory. The 
coupling between the spin-singlet order parameter t/^ and 
the spin-triplet order parameter 77 requires that we take 
spin-orbit coupling into account. The difference between 
spin-orbit coupling in Ru and Sr2Ru04 yields spin flip 
tunneling at the interface which yields a coupling be- 
tween rja and 77.^^' ^^-^ The coupling term has the form 



J's.o. = [ dS[ri*{n X J7 + c.c] (25) 

interface 

where rjs is the s-wave order parameter on the Ru-side 
and rj the two-component p-wave order parameter on 
Sr2Ru04-side, n is the interface normal vector and t is 
a coupling constant. Obviously, for n = (100) only the 
r^y-component couples. 

Further there is a coupling between the two order pa- 
rameters also away from the interface which is again due 
to spin-orbit coupling. The corresponding term in the 
free energy, derived based on symmetry arguments, has 
the form. 



^sp = J dVa[{Dxr]y - Dy-qxYns 



(26) 
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where d is again a coupling constant different on the 

two sides. This term is only active, if there is a spatial 
variation of the order parameter or a magnetic field is 
present. We find that the spatial variation of the order 
parameter along the x-axis would couple exclusively the 
77j^-component to 77^. 

This structure of order parameter coupling leads to a 
support for the nucleation of the ?7j,-component on the 
discussed [100]-interface due to the mixing with the s- 
wave order parameter intrinsic to Ru. On a qualita- 
tive level, however, the inclusion of an s-wave component 
would not modify the properties of the 3K-phase on an 
essential way. 

§6. Conclusion 

In this article we have interpreted the 3K-phase of 
Sr2Ru04 with Ru-metal inclusion as an inhomogeneous 
superconducting state located at the interface between 
the material phases. This leads to a superconducting 
phase is diifcrcnt qualitatively from the bulk phase ap- 
pearing below 1.5 K. The 3K-phase conserves time rever- 
sal symmetry. This implies an additional phase transi- 
tion where time reversal symmetry is spontaneously bro- 
ken. This second transition occurs above the onset of 
bulk superconductivity. Irregular shapes and distribu- 
tion of the Ru inclusions and the capillary effects would 
prevent a very sharp transition for the onset of the 3K- 
phase as well as the second transition. 

The origin of multiple superconducting transitions is 
analogous to that of the splitting of the superconduct- 
ing phase transition of degenerate order parameters by 
applying symmetry- lowering uniaxial stress. ^^'^^^ Indeed 
the interface represents a region of the system where the 
symmetry is effectively lower, lifting the degeneracy be- 
tween the two order parameters rjx and ry^. Therefore the 
experimental proof of the conservation of time reversal 
symmetry in the 3K-phase would be a clear confirmation 
of a having a two-component order parameter. 

We have argued that the phase structure of the pair 
wave function can lead to frustration effects in the cou- 
pling of the superconducting order parameters on differ- 
ent Ru-inclusions or different regions on the same inclu- 
sion. This is close related with the fact that this type 
of order parameter can lead to low- or zero-energy An- 
dreev bound states, since the configuration of inclusions 
forms a complex Andreev billiard system. There are var- 
ious experimental consequences due to these properties, 
some of which have been partially already investigated 
experimentally. (1) The Ru-metal inclusions provide an 
interesting way to tunnel into Sr2Ru04. Contacts via Ru 
inclusions could reveal more about the structure of the 
quasiparticle spectrum.^*) (2) The superconducting in- 
terface states of different neighboring Ru-inclusions over- 
lap and form a complex network. The study of the crit- 
ical current as a function of temperature may give an- 
other tool to investigate the unconventional nature of 
this state. In particular, the frustration effects men- 
tioned above and the spontaneous currents which oc- 
cur in the time reversal symmetry breaking phase can 
yield characteristic anomalous behavior. (3) The low- 
field magnetic response may show strong non-linear be- 



havior and cooling history-dependence in the time rever- 
sal symmetry breaking state for temperatures above Tcs , 
a phenomena which could be rather similar to the param- 
agnetic Meissner (or Wohlleben) effect in granular high- 
temperature superconductors.^^' ^^''^^^ (4) The study of 
microwave absorption in a small static magnetic field 
may show a non-monotonic field dependence as in some 
granular high-temperature superconductors.^^' 

These is only a selection of possible unusual properties 
of the 3K-phase. It is obvious that this phase and its 
properties provide a very good tool to investigate the 
superconductivity in Sr2Ru04 from a new point of view. 
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